Deconfinement transition and string tensions in SU(4) Yang-Mills Theory 
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We present results from numerical lattice calculations of SU(4) Yang-Mills theory. This work has 
two goals: to determine the order of the finite temperature deconfinement transition on an N t — 6 
lattice and to study the string tensions between static charges in the irreducible representations of 
SU(4). Motivated by Pisarski and Tytgat's argument that a second-order SU(oo) deconfinement 
transition would explain some features of the SU(3) and QCD transitions, we confirm older results 
on a coarser, Nt = 4, lattice. We see a clear two-phase coexistence signal in the order parameter, 
characteristic of a first-order transition, at S/g 2 = 10.79 on a 6 x 20 3 lattice, on which we also 
compute a latent heat of Ae « 0.6esB- Computing Polyakov loop correlation functions we calculate 
the string tension at finite temperature in the confined phase between fundamental charges, ai, 
between diquark charges, 02, and between adjoint charges 04. We find that 1 < o^/oi < 2, and our 
result for the adjoint string tension 04 is consistent with string breaking. 
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I. INTRODUCTION 

It is well-established that the dynamics of the strong 
force are described by nonabelian gauge theory with an 
internal SU(3) symmetry and matter in the fundamental 
representation. Although the full theory of QCD con- 
tains dynamical quarks, study of SU(3) pure Yang-Mills 
theory [0 provides useful physical information. For ex- 
ample, in the valence or quenched approximation where 
gauge field configurations are generated without includ- 
ing the fermion determinant in the partition function, the 
computed light hadron spectrum differs from the experi- 
mentally measured spectrum at the 10% level j^]. This is 
convenient since Monte Carlo calculations require enor- 
mous computational effort to include dynamical quark 
effects, so many studies are done in the quenched ap- 
proximation in the interest of practicality. In the present 
paper we are interested in the phase diagram of QCD in 
the temperature-quark mass plane, and so the study of 
pure Yang-Mills theory covers the m q = 00 line. 

The confinement-deconfinement transition of QCD at 
high temperature (T i=a 100—300 MeV) has been studied 
with and without dynamical quarks and depends strongly 
on the number of light quark flavors and their masses. 
Pure gauge theory is recovered in the limit of infinite 
quark masses. In this limit the order parameter of the 
deconfinement transition is the Polyakov loop 0] , 



L{x) 



1 Nt 
— Tr Y[U (x,t). 
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In the confined phase {L)% = 0, but in the deconfmcd 
phase the Polyakov loop acquires a nonzero expectation 
value, spontaneously breaking the global Z(N C ) center 
symmetry. Since a Z(3) symmetry admits a cubic term 
in the effective potential, it drives a first-order transition; 
such is not the case for N c > 3 B. 



For Nf flavors of massless quarks, the QCD Lagrangian 
has a global S\J(Nf) L ®S\J(Nf) R chiral symmetry. At 
zero temperature this symmetry is spontaneously bro- 
ken and the pions are the massless Goldstone bosons, 
but at some finite temperature T x the chiral symme- 
try is restored. Universality arguments suggest that the 
transition should be first-order for Nf > 3 ||, while a 
second-order transition for Nf — 2 is not ruled out || 
(in fact a second-order transition is supported by lattice 
studies Q). In nature the strange quark mass, m s , is 
roughly 25 times larger than the average of up and down 
quark masses, m u ,d = {jn u J rm,i)l c l. According to lattice 
calculations, the order of this "2+1" flavor phase tran- 
sition depends on the strange mass. As m s is increased 
from m Ul d, the first-order phase transition weakens into 
a crossover [|| . In Fig. |l| we reproduce the "Columbia" 
phase diagram which shows the order of the transition 
for different regions in (Tn U)C 2,rn s ) plane. 

Recently Pisarski and Tytgat [jllj argued that the 
Columbia diagram is hard to understand in light of in- 
tuitive large-iVc arguments. They point out that since 
anomaly effects are suppressed by 1/N C , the contribution 
of chiral symmetry restoration to the free energy is 0{N C ) 
while the change in the free energy due to deconfinement 
is O(Nc). So, in the large- N c limit a first-order decon- 
finement transition should be robust for any quark mass. 
Thus, if the the first-order transition of SU(3) is a general 
feature of SU(iV c ) it is hard to understand why it disap- 
pears as the quark masses increase away from zero. One 
resolution of this conflict, they propose, is that N c = 3 
is special due to the cubic term in the effective potential, 
and that the general SU(iV c ) deconfinement transition is 
second order. 

Yang-Mills theory with N c colors, SU(A C ) pure-gauge 
theory, has been a topic of exploration for lattice Monte 
Carlo study since the first days of the field . A first- 
order phase transition in the average energy was observed 
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FIG. 1: Columbia diagram showing the nature of the 2+1- 
flavor finite temperature transition for different values of m Uy d 
and m s . Regions labeled I have a first-order transition, while 
the region II has only a crossover. For m u = mj = there is 
a tri-critical point (square) above which the chiral transition 
is second-order, denoted by asterisks. The pure glue limit is 
indicated by the octagon at (oo, oo). The question mark indi- 
cates the physical (m U: d, m s ) according to calculations with 
staggered fermions |)|. (Ref. [[[o| suggests that the physical 
point lies lower, in the Region I, when Wilson fermions are 
used.) 



on symmetric lattices with volumes between 3 4 and 6 4 for 
N c = 2 - 5 [l4| H§ |l6j. Refs. |7|, f[§ gl| demon- 

strated this transition is the consequence of a lattice- 
induced critical line separating the strong and weak cou- 
pling regimes. Specifically, they added to the usual fun- 
damental single-plaquette action 







X fM,W.fl<L' 



the adjoint action 



Sa = 



6 ^ 



E 



Re Ti A P^ix), 



(2) 



(3) 



where P^ u (x) is the plaquette Ul(x)UUx + v)U v (x + 
/t)C/ M (x) and Try and Tr^ are traces in the fundamen- 
tal and adjoint representations, respectively; \x and v 
are space-time (or space-temperature) indices. Figure ^ 
shows the resulting phase diagram for this mixed action 
for SU(3) and SU(4). A first-order transition line sepa- 
rates the strong and weak coupling regimes of the funda- 
mental coupling (3. For iV c < 3 this line ends in a critical 
point before crossing the /3a = axis, but does cross this 
axis for N c > 3. Also displayed is the (3 ss transition 
line corresponding to the transition in SO(N^ — 1) gauge 
theory. 

This lattice-induced (Ad) bulk transition, while inter- 
esting, can obscure the physical (3d) finite temperature 
transition of interest here. If the two transitions are 



FIG. 2: Phase diagram of fundamental- adjoint lattice action 
showing the lattice-induced bulk transition. The transition 
line crosses the fundamental axis for N c > 3. 



nearby in parameter space, the first-order nature of the 
bulk transition would have a non-trivial effect on the 
confinement-deconfinement transition. The bulk tran- 
sition was found to occur near f3 = 10.2 (with (3 A = 
0) Jl7[. The past finite temperature studies of SU(4) 
p0|, pl|, p2| p3[ addressed this issue to varying degrees. 
For example, Ref. |2^] studied the deconfmement transi- 
tion along the line (3a = —0/2, hoping to avoid crossing 
the bulk transition line. Ref. |E3| found evidence on small 
volumes and w 1000 Monte Carlo evolution sweeps for 
deconfinement transitions for both N t = A and N t — 5 
around f3 = 10.5 and (3 = 10.6 (with fi A = 0). All these 
studies concluded that the SU(4) deconfinement transi- 
tion was first-order. However, given the alluring expla- 
nation for the Columbia diagram, we felt the time was 
ripe to revisit finite temperature SU(4) numerically. 

Another topic which we address in this work is the 
tension of confining strings which carry k £ {1,...,A C } 
units of flux. Only for N c > 3 can one find differ- 
ent strings with unequal tensions. With these finite- 
temperature calculations we find the ratio of diquark 
(ffc=2) to fundamental (<Tk=i) string tensions to be in the 
range 1 < Ok=2/vk=i < 2. As pointed out in Ref. 
these types of computations may test dualities between 
gauge theories and string theories. String tensions can 
be computed on the lattice, in broken supersymmetric 
gauge theory, and in M theory versions of QCD and su- 
persymmetric QCD. Our result for CT2/C1 indicates that 
in SU(4) Yang-Mills flux tubes attract each other as ex- 
pected from SUSY Yang-Mills and M theory (24) and 
proved in standard Yang-Mills [^5] . 

In the next Section we describe some details of our 



calculations. Section III gives the results for the decon- 
finement transition, and Section IV shows our calculation 



of the string tensions. Finally we summarize our results 
111 Section [y|. 



3 



II. COMPUTATION 



Our calculations of SU(4) Yang-Mills theory do not 
differ significantly from standard SU(3) calculations. 
We use the fundamental single-plaquette action, Equa- 
tion (||), with (3 = 2N c /g 2 = 8/g 2 . Our production code 
is a minimally modified version of the MILC code [ p6[ . 
Our algorithm for evolving the gauge fields is a mixed 
overrelaxation/heatbath procedure: in one Monte Carlo 
"sweep" we perform 10 microcanonical overrelaxation 
steps followed by one Kennedy-Pendleton [^7j heatbath 
step. Each sweep we compute the average plaquette 
and fundamental Polyakov loop. We generate at least 
1000 Monte Carlo sweeps at each (3, with 10000 to 20000 
sweeps around (3 C . An independent Metropolis code 
was written from scratch for SU(4) to check this MILC- 
derived SU(4) code. 

For those values of the coupling (3 where we want to 
calculate the string tension, we compute correlation func- 
tions of Polyakov loops in the irreducible representations 
of SU(4) 
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FIG. 3: Magnitude of the Polyakov loop vs. f3 on a 4 x 8 3 
lattice. 
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d(r) = {Li{x)L*{S + r)) s , 



(4) 



where i =4, 6, 10 and 15 and the trace in Eq. ([!]) is i- 
dimensional. As is well known, the diquark 6 and 10 rep- 
resentations are obtained by antisymmetrizing and sym- 
metrizing two fundamental 4 representations, and the 
adjoint 15 by inserting SU(4) Gell-Mann matrices. We 
use the Parisi-Petronzio-Rapuano multihit variance re- 
duction method [^] to reduce noise. Polyakov loop cor- 
relation functions are computed every tenth Monte Carlo 
sweep. We investigate autocorrelations by including only 
every n-th configuration, where n = 1, 5, and 10. We 
have a total of 2800 measurements for the calculations 
on a 6 x 16 3 lattice and 1900 measurements for the cal- 
culations on a 8 x 12 3 lattice. 



III. DECONFINEMENT TRANSITION 
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FIG. 4: Plaquette vs. /3 on 4 x 8 lattice. The normalization 
is such that ([/ p i aq ) = N c in the free theory. 



the same (3 (10.2-10.6) region as for N t = 4 (see Fig. ||), 
as is expected for a bulk transition. And with N t = 6 the 
bulk and finite-temperature phase transitions are clearly 
separated. 

The six plots shown in Fig. g show the real and imagi- 



In order to make contact with previous finite tempera- 
ture SU(4) calculations, we compute thermodynamic ob- 
servables on a 4 x 8 3 lattice for values of (3 = 8/g 2 be- 
tween 10.0 and 10.6. We find a rapid change in (\L\) 
between (3 = 10.45 and j3 = 10.5 (see Fig. |), in agree- 
ment with Refs. [^0| |2j|. Since the plaquette is also in- 
creasing in that region (Fig. Q), one might worry that the 
bulk transition is affecting the deconfinement transition. 
Therefore, we did not pursue confirming the order of the 
deconfinement transition with Nt = 4, and instead focus 
on Nt = 6 where the bulk and deconfinement transitions 
should be further separated. 

In our 6x 12 3 calculations we find that the jump in (\L\) 
is between = 10.75 and (3 = 10.80 (Fig. g). This move 
in f3 of the critical point is consistent with the conjecture 
that the deconfinement is a thermodynamic phenomenon. 
On the other hand the increase in the plaquette is over 
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FIG. 5: Magnitude of the fundamental Polyakov loop vs. (3 
on a 6 x 12 3 lattice. 
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FIG. 6: Plaquette vs. (3 on a 6x 12 lattice. The normalization 
is such that (L/pi aq ) = N c in the free theory. 
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FIG. 7: Plots at six values of /3 showing the phase and magni- 
tude of the average fundamental Polyakov loop on consecutive 
6 x 12 3 configurations. Each cross corresponds to the value 
of L4, computed on a single configuration. 



nary parts of the Polyakov loop for the last 2000 sweeps 
of calculations at the corresponding couplings. One can 
see the spontaneous breaking of the Z(4) symmetry as 
the deconfinement transition is crossed. 

One quantitative estimation of the critical coupling, /3 C , 
comes from the deconfinement fraction . Let us define 
4> £ [0, 7r/4] as the angle between arg(L4) and the nearest 



FIG. 8: Deconfinement fraction (defined in text) vs. /3 for the 
6 x 12 3 lattice. 
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FIG. 9: Deconfinement fraction vs. (3 for the 6 x 16 3 lattice. 

Z(4) symmetry axis. Given another angle 9 S [0,7r/4], 
one counts the number of configurations where 4> < 9, 
Ni n , versus the number of configurations where 4> > 8, 
N out . In the confined, Z(4)-symmetric phase, on average, 

N„, 6 



(tt/4) 



(5) 



The deconfinement fraction is the excess number of con- 
figurations which have < 9: 



m = 



tt/4 



(tt/4) 



N in + N out 



(tt/4) 



(6) 



where the factor outside the brackets normalizes the to- 
tally deconfined f(0) to one. (Note that due to statistical 
fluctuations, f(9) can be slightly negative in the confined 
phase.) The critical coupling is defined to be the value 
of (3 for which f(9) = 1/2. In Fig. || we plot the decon- 
finement fraction with 9 = 15° and 25° for the 6 x 12 3 
lattice. We find (3 C — 10.78 ±0.01, where the uncertainty 
is estimated by varying 9 between 15° and 30°. 

In order to determine the order of the phase transition, 
we increased the spatial volume to 16 3 and 20 3 . With the 
larger volumes, the critical coupling increases slightly to 
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FIG. 10: Histogram of [L 4 | at (3 C = 10.78,10.79,10.79 on 
volumes of 12 3 , 16 3 , 20 3 , respectively. 



[3 C = 10.79 as is expected (see Fig.^). The histograms of 
Polyakov loop magnitude \ (L4) \ obtained from the larger 
two lattice volumes near their respective critical points 
show two peaks, in clear contrast to the 12 3 volume. See 
Figure 10. This suggests a first-order phase transition. 
Indeed Polyakov loop evolution in simulation time, in 



Figure |Tl|, signals coexistence of the confined and decon- 
fined phases at this temperature, — 10.79. The mag- 
nitude stays with its low (confined) or high (deconfined) 
value for a relatively long period, but occasionally jumps 
very quickly from one to the other value. And when the 
magnitude is low, the argument takes random arbitrary 
values, while it is fixed to the neighborhood of one of the 
four allowed Z(4) values when the magnitude is high. 

By combining these histogram and evolution observa- 
tions, we conclude that the finite-temperature deconfin- 
ing phase transition of SU(4) Yang-Mills system is of 
first-order. It is thus desirable to compute the latent 
heat through combinations of the energy density e and 
the pressure p BO] . Specifically we compute 



a 4 (e-3p) = 
a 4 (e+p) = 



-6N c a 



dg 
da 



-(Pt+Ps 



8Nc 
9 2 



■C(g 2 )(P t -P s ). 



(7) 
(8) 



The average space-space and space-time (or space- 
temperature) plaquettes are normalized such that if g 2 = 

o,p7 = p7 = i: 



x i 

P° = i^E E ReTr^x) (10) 



x i,j:i<j 



where Q is the 4c? volume and i and j are spatial indices. 
In bare lattice perturbation theory the (3— function and 
Karsch coefficient are given, respectively, by |31j] 



da 



1F/V C 
! 3(16tt 2 ) 



0( 5 2 



(11) 
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FIG. 11: Monte Carlo evolution of the magnitude (top) and 
the argument (bottom) of the fundamental Polyakov loop for 
P = 10.79 on the 6 x20 3 volume. We take arg(L 4 ) £ 
for clarity. 
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FIG. 12: Energy density minus 3 times the pressure vs. f3 for 
the 6 x 16 3 lattice. A divergent vacuum contribution remains 
to be subtracted. 



and 

C(g 2 ) = 1 



4N r . 
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0.586844- 0.005306 ) 



(12) 



6(16tt 2 ) 

It is possible, and advisable, to use mean-field improved 
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TABLE I: Discontinuities in e — 3p and e + p at j3 = 10.79 for 
various values of r (described in text). The lattice spacing has 
been set through the critical temperature: a([3 c ) — (7VtT c ) _1 . 
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FIG. 13: Energy density plus pressure vs. f3 for the 6 x 16 3 lat- 
tice. Squares correspond to averaging over the whole data set, 
the fancy crosses at the critical beta correspond to separating 
into hot and cold phases as described in the text. 



perturbation theory or a nonperturbative calculation of 
these quantities for an accurate calculation of the energy 
and pressure p2 |; however, for the purpose of establishing 
a nonzero latent heat, bare perturbation theory suffices. 

The quantities e — 3p and e+p are plotted as functions 
of j3 on the 16 3 lattice and shown in Figures |lj and |l^, 
respectively. The fancy crosses in the latter figure cor- 
respond to separating the configurations at (3 C = 10.79 
into hot and cold phases. Note that e — 3p plotted in 
Fig. [12] contains a divergent vacuum contribution which 
may be subtracted after a zero temperature simulation is 
performed; however, such subtraction is not necessary in 
order to compute the latent heat from a discontinuity in 
e— 3p at the critical coupling (3 C . The separation of phases 
at C was made on the basis of whether \L±\ was greater 
or lesser than some value r. Based on the histograms in 
Fig. |h] we varied r from 0.08 to 0.14. Table § lists the 
values for A(e — 3p) and A(e + p) obtained for different 
r on both the 16 3 and 20 3 volumes. The variation as a 
function of r is within the statistical errors. 

Thus, we observe a latent heat which is many standard 
deviations greater than zero. We also see that A(e— 3p) ^ 
A(e+p) which implies a discontinuous change in pressure 
across the transition. If, for example, we take the 6 x 20 3 
data with r — 0.10, we find (with statistical errors only) 



Ae = 5.7(3)T C 4 
Ap = -0.45(13)T C 4 . 



(13) 
(14) 



A nonzero Ap was also seen in early studies of SU(3) 
p3| , |34| and disappeared when going from the perturba- 
tive estimates for ( |ll[ ) and (JlJ) to nonperturbative cal- 
culations p2"| |. 

A thorough calculation of the latent heat in the SU(4) 
deconfinement transition requires a full study of the lat- 
tice spacing dependence as well as nonperturbative deter- 
mination of the /3-function and Karsch coefficient. How- 



ever, even the exploratory study here makes clear the 
latent heat is nonzero and further establishes the first- 
order nature of the phase transition. 

We can compare the latent heat for SU(4) to that for 
SU(3) by normalizing by the energy density for an ideal 
gluon gas. If we take the latent heat to be Ae/T c = 6.0 ± 
1.5 and divide by the Stefan-Boltzmann energy density, 



we find 



e SB (T) = 



Ae 

esB(T c ) 



(ni - iy< 

15 



rp4 



0.60 ±0.15. 



(15) 



(16) 



Our result should be compared against the N t = 6 SU(3) 
latent heat obtained using a perturbative f3— function: 
Ae/esB = 0.454(11) §5). A state-of-the-art SU(3) calcu- 



lation, which used an improved action and a nonpertur- 
bative ^-function, gave Ae/e SB = 0.266(17) ff|. Fur- 
ther work is required to see if the effect of going from 
a perturbative to nonperturbative (3— function is as dra- 
matic for SU(4) as for SU(3). 



IV. STRING TENSIONS 

We use two different lattices, 6 x 16 3 and 8 x 12 3 , for 
studying string tensions. For the former, we choose the 
coupling values of j3= 10.65 and 10.70, safely away from 
both the bulk and the deconfining phase transitions, and 
in the confining phase (see Figure Q). Polyakov loop cor- 
relations (see Eq. (ji|)) for the fundamental (4, k = 1, 
top) and anti-symmetric diquark (6, k = 2, bottom) rep- 
resentations are shown in Figures [li] and [l^. A clear 
difference in the rates of exponential decay is observed 
between C4 and C 6 - Using a correlated, jackknifed fit to 



the form |36 



a k 



exp[-V k (r)N t ] 



■exp[-V k (N.-r)N t ], (17) 
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FIG. 14: Polyakov loop correlation function in 4 (top) and 6 
(bottom) representations on a 6 x 16 3 lattice at /3 = 10.65. 
The symbols are the data points and the solid lines are fits in 
the range 4 < r < 12. 
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FIG. 15: Polyakov loop correlation function in 4 (top) and 6 
(bottom) representations on a 6 x 16 3 lattice at f3 = 10.70. 
The symbols are the data points and the solid lines are fits in 
the range 5 < r < 11. 



with N t = 6, N s = 16 and 



ments. However, we can infer from our data that 



Vk(r) = (r k r - 



3N? 



(18) 



we obtain string tensions, a% and 02, and their ratio, 
tabulated in Table ||. The analysis of the correlation 
functions is done using every measurement, every fifth 
measurement, and every tenth measurement in order to 
estimate correlations between successive measurements 
(each separated by 10 Monte Carlo steps, see Sec. ||). 
The increase in the statistical error with the number 
of skipped configurations, iV s kip, indicates a significant 
auto-correlation. Unfortunately, it appears that several 
hundred configurations are necessary in order to obtain 
a precise fit, so we cannot drop too many of the measure- 



TABLE II: String tensions (in lattice units) between static 
fundamental (k = 1) and diquark (k — 2) charges, and their 
ratio, using every (JV s ki p + l)th configuration in the analysis. 
The quoted uncertainties are statistical. 






iVskip 


(Tl (T2 (72 l<J\ 







0.098(2) 0.138(14) 1.45(15) 


10.65 
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0.092(4) 0.137(30) 1.67(36) 




9 


0.101(7) 0.164(56) 1.77(67) 







0.076(2) 0.118(13) 1.59(13) 


10.70 


4 


0.080(4) 0.154(40) 1.91(49) 
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0.084(6) 0.163(53) 2.03(68) 
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CI 



> 1 



(19) 



by roughly 2 standard deviations. Note that both a\ 
and (T2 decrease as j3 — > /3 C (i.e. as T increases). Since 
the lattice spacing decreases as T increases, the fit range 
for f3 = 10.70 does not include the r = 4 and r = 12 
data (see captions of Figs. [l4] and |l5|). Our numerical 
accuracy is good enough to conclude there are two dif- 
ferent strings, one between the fundamental charges car- 
rying one unit of flux, and another, stronger, between 
the diquark charges carrying two units of flux. It is not 
yet good enough, however, to distinguish among various 
predictions for this ratio summarized by Strassler pl[ . 
However, this establishes numerically the expectation for 
(Ti ^ <t 2 in SU(4) Yang-Mills theory, just as Ref. |37j 
showed (Ti = o"2 in SU(3) Yang-Mills theory. 
A string model |3H] predicts that 



y/cn(T = 0) \<d-2) 



= 0.69, 



(20) 



which is quite close for SU(3) H|. We have not computed 
the zero temperature string tension, but only the string 
tension roughly near T c , to find 



vVi(T«T c ) 



= 0.60. 



(21) 



The extent which the lattice scale changes between (3 — 
10.70 and (3 C = 10.79 is main uncertainty above. Of 
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FIG. 16: Polyakov loop correlation function in 4 (top), 6 
(middle), and 15 (bottom) representations on a 12 x 8 lattice 
at p = 10.85 



course a zero temperature study is necessary before one 
can assess the agreement with Eq. (20). 

On this lattice of 6 x 16 3 , which is coarser and larger 
of the two, no signal was obtained for either the sym- 
metric diquark (10) or adjoint (15) representations. In 
contrast, with the finer lattice spacing (at (3 = 10.85) 
on the smaller 8 x 12 3 lattice, flattening of the adjoint 
correlation is observed (see Figure |l6|) . This suggests the 
breaking of confining string for the adjoint representation 
at a rather short distance of 3 lattice spacings. It gives 
us confidence that the correlations on the 6 x 16 3 lattice 
should be dominated by the non-perturbative strings for 
ranges longer than at least 3 lattice units. Notice also 
that while string breaking is an expected behavior for 
the adjoint representations in general |||, such an ab- 
sence of string is yet to be observed in SU(3) Yang-Mills 
theory which employs much finer and larger lattices than 
the present work. 



V. CONCLUSIONS 

We have revisited the confinement-deconfinement 
transition of SU(4) Yang-Mills theory through Monte 



Carlo lattice calculation. One problem with the earlier 
results is that the deconfinement transition with N t = 4 
is very close in coupling constant space to a known bulk 
transition, so that its finite-temperature nature or its 
order is not clear. We have shown that by decreasing 
the lattice spacing by 2/3, the deconfinement transition 
moves upward in the coupling and proves itself as a finite- 
temperature transition, and it becomes well-separated 
from the bulk transition which does not move. Neverthe- 
less, we observe a clear signal for coexistence of confined 
and deconfined phases at this deconfinement transition. 
Therefore, we confirm that the deconfinement transition 
of SU(4) Yang-Mills theory is first-order. Additionally 
a first calculation of the latent heat of the SU(4) de- 
confinement transition has been presented here, giving 
Ae ~ 6T C , or Ae/egB ~ 0.6. Using improved techniques, 
the SU(3) latent heat is Ae/e SB = 0.266(17) |||, and it 
will be interesting to see how the latent heat depends on 
N c . 

Our calculations of the string tensions are a first study 
in lattice SU(4) and should be improved to meet the cur- 
rent state-of-the-art which exists for SU(3). Even so, we 
observe a ratio for 4 and 6 dimensional string tensions 
which is between 1 and 2. It also appears that the ad- 
joint string breaks at a short distance. We hope this work 
shows that it is interesting and feasible to study ratios of 
string tensions for N c > 3 lattice simulations. 
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